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CHAPTER I 


INTRODUCTION 

Wind is an important consideration in the analysis of airplane 
flight in the atmospheric boundary layer, both because of short-scale 
gusts or turbulence and because of large-scale variations of the mean 
wind. In the planetary boundary layer the mean wind decays toward the 
ground and has considerable horizontal variations due to irregulari- 
ties in terrain. Thus, both spatial and temporal variations occur in 
near surface winds encountered along ascending and descending flight 
paths , 

Previous analyses of airplane motion that have been carried out 
[1 , 2, 3] consider, in general, only constant winds and thus neglect 
the effects of wind shear. This report, however, investigates the 
influence of variable mean wind fields and discrete gusts on the 
dynamics of aircraft during terminal flight operations. Mathematical 
models of the winds are introduced into the equations of aircraft 
motion, both with fixed and automatic controls; computer solutions of 
the resulting motion are carried out. 

As an aircraft descends on its glide slope, a sudden change in 
horizontal wind or vertical wind, or both, will instantaneously affect 
the velocity of the aircraft relative to the air mass. If the shear 
is such that the relative velocity of the aircraft increases, the lift 
force will increase and the aircraft will tend to rise above the glide 
slope. If the shear causes a sudden decrease in the relative 


velocity, the aircraft will respond by falling below the glide slope, 
and a hazardous condition may result. 

Several reports have been published which link short and long 
touchdowns to a sudden wind shear occurrence during final approach 
[4, 5, 6, 7], Recent accident reports also have found wind shear to 
be at least a contributing cause for several accidents [8, 9]. In 
addition, it is believed that wind has been responsible for many other 
accidents, though it remained undetected at the time [10, 11]. The 
problem of quantitatively defining the effect of shear of given magni- 
tude on an aircraft during descent has not been completely resolved. 
Noteworthy studies that have investigated wind shear and/or turbulence 
during landing include References 3, 6, 7, and 10 through 14. 

Although a very complete development of the system of equations 
governing airplane motion is available [1, 2], most analyses reported 
to date reduce the equations to those for a constant wind or employ a 
linearized model which requires the assumption of a uniform wind field 
and is not applicable for non-uniform winds. 

Ramifications of the airplane motion due to the effects of 
temporally and spatially varying mean winds are studied in this 
report. Analyses of flight paths through changing mean wind fields 
reported in the literature are primarily two-dimensional and deal only 
with vertically varying horizontal winds (i.e., having a component 
parallel to the flat earth only). 

Etkin [1] has a very complete development of the general 
equations of unsteady motion. Normally, the wind components of 
velocity are not included in the equations since it is assumed that 
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no wind is present. Luers and Reeves [7] developed a system of 
equations in two dimensions which incorporate only horizontal wind. 
Later in this report, a general form for the two-dimensional equations 
of motion is developed. This accounts for both vertical and horizontal 
mean wind components with both time and spatial variations. 

Using this later set of equations, analyses both with fixed and 
automatic feedback control were carried out. In the fixed control 
simulation, the aircraft is trimmed at an altitude of 91 m and on a 
glide slope of -2.7 degrees. The corresponding throttle setting and 
elevator angle setting are then fixed for the remainder of the 
landing. These fixed control landing simulations were carried out for 
several different wind fields and the deviations in the glide slope 
and touchdown points are compared. In the cases with high wind shear, 
the deviations are very large and in some cases the aircraft 
trajectories with fixed controls are not realistic. 

To overcome this difficulty an automatic control system was 
developed for the same two-dimensional system of equations. Every 
phase of the flight of an aircraft can be regarded as the accomplish- 
ment of a set task, i.e., flight on a specified trajectory. That 
trajectory may simply be a straight horizontal line traversed at 
constant speed or it may be a turn, a transition from one symmetric 
flight path to another. All of these situations are characterized by 
two common features, namely, the presence of a desired state and the 
departures from it, designated as errors. These errors are, of 
course, a consequence of the unsteady nature of the environment. 
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The correction of errors requires a method of measuring the error 
or the desired state. Some of the state information needed (i.e., air- 
speed, altitude, rate-of-climb, heading, etc.) is measured by standard 
flight instrumentation. This information is not generally sufficient, 
however, when both guidance and altitude stabilization are considered. 
For this case, the state information needed may include [1, 16] 
position and velocity vectors relative to a suitable reference frame, 
vehicle altitude, aerodynamic angles, etc. A wide variety of devices 
are used to measure these and other variables, and range from pitot- 
static tubes to sophisticated inertial-guidance platforms. Gyro- 
scopes, accelerometers, magnetic and gyrocompasses, angle-of-attack 
and sideslip vanes, and other devices, all find applications as 
sensors. The most common form of output is an electrical signal, but 
fluidic devices [17] are receiving increased attention. 

In this study we assume that the desired variable can be measured 
independently and linearly, which is of course, an idealization. 

Since every sensing device, together with its associated transducer 
and amplifier, is itself a dynamic system with characteristic 
frequency response, noise, nonlinearity and cross-coupling, these 
attributes cannot be ignored in the final design of real systems, 
although one can usefully do so in preliminary work [1]. 

In the automatic control simulation, the aircraft is trimmed 
initially at an altitude of 91 m on a straight horizontal line 
trajectory and is automatically controlled by actuation of thrust and 
elevator angle. In this first phase the aircraft remains in an 
altitude hold mode until it intersects the Instrument Landing System 
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(ILS) guidance beam, It then switches to the glide-slope capture mode 
which actuates the thrust and elevator controls so as to capture the 
-2.7 degree glide path specified by the ILS guidance beam. As soon as 
the specified glide path is captured, the third phase, glide-slope 
tracking mode, becomes effective. In this mode the controls are 
actuated such that the aircraft remains on the glide path. At an 
altitude of approximately 18 m, flare initiation altitude along with 
other necessary parameters are calculated to begin the flare mode. 

The flare mode is switched on as soon as the aircraft reaches the 
flare initiation altitude. The aircraft remains in this mode until 
the final touchdown. 

For this investigation, the scope of the automatic landing 
problem was restricted in two ways. First, the aircraft simulation 
equations are restricted to three degrees of freedom by considering 
the longitudinal axis only. This restriction is reasonable in the 
light of the accident statistics compiled in References [8, 9], 
which conclude that accidents due to longitudinal errors are fatal 
more often than accidents due to lateral errors. Second, the system 
guidance information was assumed to come from error-free sensors and 
an error-free ILS beam. This is beneficial to maintain simplicity of 
the automatic control subroutines, since the objective of this study 
is the effects of wind shear, and not a study of ILS system errors. 
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CHAPTER II 


AIRCRAFT LANDING MODEL 


1. Equations of Motion 

The two-dimensional model for aircraft motion presented in this 
section follows the general form developed by Frost [12]. It accounts 
for both vertical and horizontal mean wind components having both time 
and spatial variations. 

The aircraft trajectory model employed in this study was derived 
based on the following assumptions: 

a) The earth is flat and non-rotating. 

b) The acceleration of gravity, g, is constant (9.8 m/sec ). 

c) Air density is constant (1.23 kg/m ). 

d) The airframe is a rigid body. 

e) The aircraft is constrained to motion in the vertical plane. 

f) The aircraft has a symmetry plane (the x-z plane). 

g) The mass of the aircraft is constant. 

h) Initial flight conditions are for steady-state flight. 

Figure 1 illustrates the forces acting on the aircraft. These 

i nclude: 

F-p thrust of the engines 

L lift 

5 drag 

W wind velocity 

mg gravitational force. 
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Figure 1 Relationship between the various forces acting on an 
aircraft [12] 
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The figure shows the orientation of the forces with respect to 

->■ 

the velocity relative to the earth (V), the velocity relative to the 
air mass (V), and the fuselage reference line (FRL) of the aircraft. 

a 

The x-axis in Figure 1 is parallel to the surface of the earth and the 
z-axis is perpendicular to the surface of the earth (positive down- 
ward) . 

*>• 

From a direct force balance along the direction V and along the 

4 - # 

direction perpendicular to V, respectively, it follows from Figure 1 
that 

mV =-LsinS-D cos 5 - mg sin y + Fj cos (6y + o) (1) 

and 

mVy = L cos y - D sin 6 - mg cos S + Fy sin (5y + a) . (2) 


The aerodynamic forces and the thrust from the engines exert a 
pitching moment on the aircraft. The equation describing the momentum 
balance about y is 


q = 


dt 2 


f t l t 

l yy 


M 


yy 


(3) 


where the dot refers to the derivative with respect to time, and 
g is the magnitude of the acceleration of gravity, 

V is the magnitude of the velocity relative to the earth, 

y is the angle between V and the x-axis (the flight path 

angle) , 

Fy is the magnitude of the thrust vector, 
m is the aircraft mass, 
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6 


T 


a 



is the angle between the thrust vector and the fuselage 
reference line (FRL), 
is the angle between V and the FRL, 
is the angle between and V, 

a 

is the time derivative of the pitching rate, q, 
is the effective moment arm of the thrust vector, 
is the pitching moment, and 

is the moment of inertia about the symmetry plane of the 
ai rcraft . 


By considering a different coordinate system in which the x-axis 
is along the vector V a , called "wind" frame of reference by Etkin [1], 

a 

similar force equations can be developed by summing up the forces 
parallel and perpendicular to V . These are 

m(V+W)+mqW “Ft - D - mg sin y' (4) 

a x w z T 

w w x if 

w 

and 

mW - mq (V_ + W v ) = F t - L + mg cos y' . (5) 

Z W a A, _ I _ 

w w z 

w 


It is convenient to express these in terms of the wind 
components relative to an earth fixed coordinate system, since most 
wind correlations from the meteorological literature are expressed in 
such coordinates. 

W v = W cos y’ - W sin y' , (6) 

V x 

W z = W x sin y* + VI z cos y' • (7) 
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Taking the time derivative of W , we get 

x w 


• « 


W v = W cos y' - W sin y 1 - W siny 


x w x 


, dy!_ 
dt 


W Z C0S Y' 


dy 1 

Then, since 9 W = » 


• * 


W x = W cos y 1 - W z sin y' - W q . 
w w 


Also, since 


Fj = F-j- cos( 6 T + a') , and 

x ‘ ! 

w 


Fj = F-j- sin(Sy + a 1 ) , 


Equation (4) becomes 

mV a = Fy cos(5y + a') - D - mg sin y 1 - m(W x cos Y 1 “ sin y‘ 

From Equation (7), taking the time derivative of W , we get 

w 

W z = W x sin y' + W z cos y' + W x cos y' - W z sin y 1 , 
w 

and Equation (5) becomes 

- mV a q w = - Fy s i n ( 6 y + a') - L + mg cos y' 

• • 

- m(W sin y* + W cos y') . 

A L- 

The moment equation remains the same as Equation (3). The 
governing force equations in "wind" frame of reference are thus 


) 


(8) 


(9) 


• ( 10 ) 


( 11 ) 


( 12 ) 
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mV a = Fy cos C^y + ol { ) D - mg sin y' 

- m(W cos y 1 - W sin y') , 

A Z 

mV a y' = Fy sin(Sy + a*) + L - mg cos y* 
+ m(W sin y 1 + W cos y * ) , 

A Z 


f t l t 

% = — + 

yy 


H 

x * 

A yy 


(13) 


(14) 

(15) 


where W is the horizontal wind speed, W is the vertical wind speed, 

A iL 

and a 1 (the angle of attack) is the angle between V and the FRL, 

a 

2. Incorporation of Wind Shear 

The wind is seen to enter the equations in the form of a gradient 
or wind shear W and W . The expanded form of these equations is: 

j\ Z 


W x = 
or 


5c 

at 


aw 


x dX 


aw 


x dZ 


ax dt dt 


w.. = 


aw aw aw 

x + V [ cos y - sin y ] , 


at 


(16) 


and, similarly, 

aw „ aw aw 

W Z = ^r +v[cos sin Y 1T ] • (17) 

Thus, both spatial variations and temporal variations in atmospheric 
motion influence the equations in the wind coordinate system. 


11 


3W 3W 

X z 

Generally, care is needed in evaluating and since the 

3 L oL 

wind speed is normally expressed in terms of altitude measured upward 
from the surface of the earth, whereas in aerodynamic coordinates, Z 
is measured downward. 

Additional kinematic relationships necessary to solve for the 
aircraft motion are as follows: 

The relative velocity as a function of inertial and wind velocity 
is 

9 . o 1/2 

V = [ (X - w ) 2 + (Z - W ) 2 ] , (18) 

a X l 

and, in turn, 

2 

V = W cos y - W sin y + [ (Ia/ sin y - W cos y) 

x z z a 

999 1/2 

+ V 2 - (W Z + W Z ) ] . (19) 

a x z 


->• 

The angle between V and is given by 

W sin y + W cos y 

sin 6= • 

v a 

Other angular relationships are 
a^Q-y-S^O-y 1 , 
a = 0 - y . 

The derivative of a' is 
• . X *i 

a =0-Y=q-y , 


( 20 ) 


( 21 ) 

(22) 
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where y‘ is given by Equation (14), hence, 



F-j. sin(6 T + a 1 ) 
mV a 


mV a 


, g COS Y 1 ,1 r 1*1 • , , , T 

+ * — r: — + tt- [W Sin y + W COS y' ] 


( 23 ) 


Also required for solution of the preceding equations are the 
aerodynamic coefficients 

C L = C L (a' , 6 e , V a , q, a 1 ) , 

C D ~ (a', 6^, , q, a 1 , C^) , 

C m = C m <“'• V V «>• »') > < 24 > 

where is the elevator deflection angle. As indicated above, the 
aerodynamic coefficients are functions of a number of variables. The 
expressions for C^, C^, and C , along with the stability derivative 
data and aircraft physical data are given in the Appendix. 

The equations of motion discussed in this chapter can be solved 
for the flight of an aircraft flying through spatially and temporally 
varying two-dimensional wind fields. In this study we have used three 
different wind shear models, 

1) atmospheric flow over homogeneous terrain, 

2 ) atmospheric flow over buildings, and 

3) atmospheric flow associated with thunderstorm gust fronts. 

The initial conditions used in this simulation are for that of a 

pitch stabilized aircraft, given by 
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Altitude = 91 m 

V_ = 70 m/sec 

a 

V a = 0 

Y 1 =0 

q = 0 

q =0 


a 


0 . 


Under these conditions the initial values of thrust, elevator angle 
and angle of attack were calculated from Equations (13), (14) and (15). 
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CHAPTER III 


MATHEMATICAL MODELS FOR VARIABLE WIND FIELDS 


1 . Atmospheric Flow Over a Homogeneous Terrain 

The mean velocity in the region of the atmosphere near the ground 
is described by a logarithmic function of altitude. The surface 
roughness characteristic of most natural terrains is generally 
described in terms of a vertical scale, Z Q . For a neutral atmosphere, 
experimental evidence [18, 19] confirms that the mean wind velocity in 
the region near the ground can be described by a logarithmic wind pro- 
file (Figure 2). The logarithmic wind profile is thus [20], given as 
a function of altitude Z, 


u* 

W X W = — in 


Z + Z 




(25) 


where Z Q represents the surface roughness, and k is von Karman's 
constant, u* is the friction velocity given by u* » where t q 

is the surface shear stress, and p is the air density. Observed 
wind profiles up to 150 m, over reasonably level and uniformly rough 
terrain, with neutrally stable conditions, obey this law reasonably 
well [20]. 


2. Atmospheric Flow Over Buildings 

Since the wind shear models for the flow over a block building 
and a step are completely described by Sheih, et al . [22] and Bitte 
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Figure 2 Logarithmic wind profile 
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and Frost [23], respectively, only a cursory description of the models 
is presented here. 

The distorted shear flows approaching and passing over a building 
can be divided into a displacement zone, an upstream bubble or down- 
wash zone, and a wake zone which includes the rear separation bubble 
or cavity zone (see Figure 3). The effect of shear in the approaching 
flow creates a downwash on the front face and a swirling flow in the 
wake or cavity zone. Undisturbed, neutrally stable atmospheric wind 
perpendicular to the axis of the building is assumed far upstream and 
far downstream of the obstacle (see Figures 4 and 5). The atmospheric 
wind field is analyzed by using the Navier-Stokes equations with a 
two-equation model, one for the turbulence kinetic energy and the 
other for turbulence length scale. In this approach, the partial 
differential equations for the vorticity, stream function, turbulence 
kinetic energy, and turbulence length scale are solved by a finite- 
difference technique. 

Both vorticity contours (Figure 6) and streamline patterns 
(Figure 7) confirm the experimental evidence of a small downwash zone 
near the front lower corner and a large recirculation zone behind the 
obstruction. Figures 8 and 9 show the computed velocity profiles at 
selected X-stations in the region close to the wall. The flow is 
decelerated as it approaches the obstruction and is accelerated as it 
passes over the obstruction. In the region above the recirculation 
zone, the flow is accelerated because of the displacement of the flow. 
The flow re-attaches near X = 12. 3H and the logarithmic boundary 
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Upstream separation Rear separation 

bubble or down wash bubble or cavity 


zone 



Figure 3 Definition of flow zones near a sharp-edged block [22] 
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X = -10.0 H Wall Boundary X » 20.75 H 


Figure 4 Description of flow region considered for a block building [22] 


Z = 9 . 0 H Upper Boundary 



Figure 5 Description of flow region considered for a long, wide building [23] 


Height Z/H 



3 


4.03 
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8.03 


Horizont. 


Figure 7 Stream! 



X/H 

(H = 3.2 m) 


Figure 8 Velocity profiles over an obstruction on the surface [22] 
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Horizontal Distance X/H 


Figure 9 Velocity profiles over a step geometry long wide building [23] 


layer begins to re-establish downstream. These results seem to agree 
very well with the limited experimental data available. 

The velocity distributions of the atmospheric flow around 
buildings are especially important in the design and operational pro- 
cedures for helicopters and V/STOL aircraft operating in large 
metropolitan areas, 

3^ Atmospheric Flow Associated with Thunderstorm Gust Front s 

Gusty winds are undoubtedly the most hazardous for an aircraft to 
negotiate. One of the most common causes of significant wind shear is 
the gust front associated with thunderstorms. The thunderstorm gust 
front is believed responsible for several accidents [5, 22]. The 
severe wind shear accompanying thunderstorms is generated by a 
vigorous rain-cooled downdraft, which spreads out horizontally from 
the storm cell as it approaches the ground. The cold outflow is led 
by a strong, gusty wind which often occurs as much as 16 km ahead of 
the storm, called the gust front. 

Mathematical schemes for computing wind fields associated with 
thunderstorm gust fronts are still in the formative stages. After 
extensive study of gust front characteristics and the available gust 
front data, Fichtl and Camp [26] have presented a mathematical model 
which describes updrafts and downdrafts associated with gust fronts 
along a given approach path. This model incorporates both scaled 
vertical wind speeds along a -2.7 degree glide slope from the gust 
front data of Goff [27], and the vertical wind speeds reconstructed 
from the digital flight data record of Eastern 66 [8]. The sequence 
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of vertical wind speeds encountered by an aircraft during landing is 
given by the following: 


Major downdraft: 

X - X 

W z = - P 1 A sinU-^-T) ; Z 0 > Z > Z r 
Major updraft: 

0 - 2q 0 )(X - X r ) 3 + (1 -3q2)(X-X r ) 2 + (2q 0 -3q2)(X-X r ) 

W = A ~ — — 


Minor downdraft: 

X - 1 - X 

W z = - P 2 A sin (it ) ; (Z p - L) > Z > (Z p - (1 + q 2 ) L 1 ) 

Minor updraft: 

X - 1 - x 

W z = P 2 A sin (tt ) ; [ Z p - (1 + q 2 ) L'] > Z > [ Z p - (1 + 2q 2 ) L'] 

where 



The various quantities in the above equations are defined as follows: 

W z - thunderstorm cold air downdraft, 

Z = vertical coordinate, 

Zg = altitude of the top of the major downdraft, 
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Z r -altitude of the top of the major updraft, 

Zg = altitude of the bottom of the minor downdraft, 

A = amplitude of major vertical velocity updraft, 

L' = vertical extent of major vertical velocity updraft relative 
to the flight path, 

P-j = ratio of major downdraft to major updraft velocities, 

P 2 = ratio of minor downdraft or minor updraft to major updraft 
velocity , 

% = < Z r - Z J' L '’ 

<*1 “ (Z r - z 0 )/L', 

q 2 = (Z r - Z B - L ' )/L 1 . 

The values for the cold air outflow parameters used in this simu- 
lation study are that of typical vertical wind speeds derived for data 
provided by NOAA/NSSL [27]: 

L' = 91 m 
Z R = 152 m 

A = 15.0 

P 1 = 1.2 

P 2 = 0.35 

q Q = -0.36 

^1 = 2.0 

^2 = 2.3. 

In this wind shear model, the vertical winds due to the thunder- 
storm gust front described above are superimposed on a stable 
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atmospheric boundary layer. The mean wind velocity under stable con 
ditions is given by 

u * { , Z + 1 o, 1 

WJz) = — In ( -=-2. ) + 5.2 f 

x K [ L o L 

where L is the Monin-Obukhov stability length. 
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CHAPTER IV 


AUTOMATIC CONTROL SYSTEM 

The two principal quantities that need to be controlled in 
symmetric flight are speed and flight path angle, i.e., the vehicle 
velocity. To achieve this, control forces are needed both parallel 
and perpendicular to the flight path. Parallel force is provided by 
thrust or drag control, and perpendicular force by lift control 
achieved via elevator deflection. It is evident from simple physical 
reasoning (or from the equations of motion) that the main initial 
response to opening the throttle (increasing the thrust) is a forward 
acceleration, i.e., speed control. The main initial response to 
elevator deflection is a rotation in pitch, with subsequent changes 
in angle of attack and lift, and hence, development of y, a rate of 
change of flight path angle. When the transients that follow such 
control actions have died away, a new steady state is achieved. ■ 

In this section the longitudinal automatic landing system will be 
described, and some of the design considerations will be given. 

It has been shown in Reference [24] that turbulence causes larger 
deviations from the desired flight path than the errors in ILS 
guidance. This study, therefore, concentrates on the effect of wind 
shear on safe automatic landings. The system guidance information is 
assumed to come from an error-free ILS beam and altimeter. Reference 
[3] has an excellent discussion on longitudinal automatic landing and 
aircraft control laws. The overall control system can be represented 
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in the form of a block diagram as shown in Figure 10, The flight 
control laws are segmented into control modes for different portions 
of the approach and landing, 

In the various control modes the linear filters are described by 
transfer functions. The filtering function, however, is actually per- 
formed in the digital computer by solving difference equations. 

Simulating a linear system with the appropriate difference 
equation is more efficient than solving the differential equations 
directly by numerical integration. Numerical integration would be 
rather a lengthy process and may be unstable for large sampling 
intervals [3]. The difference equation for one of the linear filters 
used in this simulation can be derived as follows: 


r ^\ = Y ( s ) - Ks 

G(s) " X(s) ~ s + a * 


The response of the filter to a unit step input is 


Y(s) = or y ( t ) = K e“ at 6(t) . 

s + a s s+a J 


Taking the z-transform, 


y(z) 



Factoring the z-transform of a unit step function from y(z), 


y(z) = K 


z-1 


z - e 


-aT 


z 

z-1 
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Figure 10 Overall block diagram of control system 





Since the response to a general input is wanted, replace the step 
input by a general input x(z). The actual input is approximated by a 
linear combination of unit step functions and is equal to the general 
input at the sampling points, 

y(z) = K - - x(z) . 
z - e~ a 

Cross-mul ti plyi ng , 

(z - e~ aT ) y(z) = K(z - 1) x(z) , 

Dividing throughout by z, 

(1 - ^y-)y(z) - K ( 1 - 1) x(z) . 

Transforming to the discrete time domain, 

y n - ^ Vl = K(x n ' x n-l } ’ 

and rearranging, the final form of the difference equation is obtained, 

' e ' aT Vl + K(x n ' Vl } * 

At the sampling points the difference equation is the exact 
solution for the response of the equivalent analog system. The 
difference equations for the various filters were compared to those 
obtained using numerical integration (fourth-order Adams-Bashford) ; it 
was concluded that the computation efficiency of the difference 
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equations is better than that of the numerical integration tech- 
nique [3]. 

1 . Mode Selector 

The mode selector automatically selects the proper control mode 
in sequence (i.e., the altitude hold, glide slope capture, glide-slope 
tracking, and flare mode) according to predetermined criteria defined 
by the desired flight path. 

In accordance with conventional practice [25], the control modes 
operate on the velocity and pitch stabilized aircraft and therefore 
operate with only two command variables: speed command, V , and 

pitch angle command, e . 

The mode selector is best described by considering a landing 
approach (see Figure 11) and the flow chart (Figure 12). The aircraft 
approaches the ILS glide slope at a constant altitude of 91 m until 
the aircraft intercepts the ILS beam. During this time the mode 
selector maintains Mode 1, i.e., the altitude hold mode. As soon as 
the aircraft penetrates the ILS beam, the mode selector compares the 
aircraft position to the point of intersection of the horizontal 
flight path and the ILS glide slope; when the aircraft reaches that 
point it switches to Mode 2, the capture mode. The capture mode is 
timed and the flight path angle is compared with the desired glide 
slope; as soon as the desired glide slope is reached the mode selector 

switches to Mode 3, the glide-slope tracking mode. Glide-slope 
tracking proceeds to a preselected altitude, at which point the sink 
rate and velocity of the aircraft are used to calculate the flare 
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Figure 11 Automatic landing geometry using ILS 
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initiation altitude (h ste p) and other initiation parameters described 
in the flare-mode selection. The mode selector then switches to 
Mode 4, the flare mode, when h stGp is reached. 

Notice that no automatic go-around mode is provided; the simu- 
lated aircraft is forced to land so that the conditions can be found 
that result in unsatisfactory landings. 


2. Altitude Hold Mode 

A simple hold mode incorporated in the system keeps the aircraft 
flying at a constant altitude. The digital control was modeled after 
the representative analog system shown in Figure 13. The system con- 
sists of a differencing circuit for calculating the altitude error, 
Ah, followed by a low pass filter, a gain, and a low gain integrator. 
For comparison, digital equivalent equations are as follows: 


e = c, e r + c 9 e r + c 3 Ah. + c Ah,, , 
C 1 c^ 2 Cj _ 2 3 j y j 1 


Ah j " h j - h ref ’ 


9 j-l + c 2 Ah j-1 


5 


These equations are solved once each computation cycle. The 
constants, Cp c 2 and c^, are given by 
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Figure 13 Altitude hold mode 
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c-j = e 


a i T 


c 2 = a, ^ ‘ c l^ ’ 


C3 ( ^2 T - 1 ) , 


where T is the sampling interval. 


3. Capture Mode 

The capture mode (Figure 14) provides for a smooth rotation from 
level flight to the glide-slope angle. In this mode a step pitch 
angle command, A9 p5 is applied to rotate the airplane. The magnitude 
of the step is based upon thq glide slope angle of the beam to be 
captured. In addition, an inertial vertical velocity error signal is 
generated to increase the sink rate for a given glide-slope angle. 

The error signal, 0 , is then integrated and filtered to produce the 
pitch angle command. The integrator provides an error signal pro- 
portional to the altitude error. Since the sink rate reference value 
is the proper sink rate of the aircraft on the glide slope, the 
resulting altitude reference is a parabolic curve that smoothly 
intersects the glide slope. 

The difference equations are as follows: 




C M 


c 2 9 c + c 3 + c 4 . 

6 C j-2 J e j-l 4 e j-2 


+ A0. 


38 




0 

C 

Figure 14 Glide slope capture mode 
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where 


c 


1 


1 + c 


-T/t 




C 3 = k^ k^ [ T + t-j (c£ “1)1 s 
c 4 = k 3 k 4 [ T " c 2 ( T + T ) ] • 


4. Glide-Slope Tracking Mode 

After glide-slope capture the aircraft remains in the glide-slope 
tracking mode (Figure 15) until flare. The glide-slope error angle, 
e, is passed through a low-pass filter, a gain, and then a low-gain 
integrator. In addition, two differencing circuits are used which 
estimate the approximate altitude error for the next step and give 
pitch error signal, 0. This extension is a pitch altitude command, 

© c , proportional to h - h r and also h-h c , where h r is the reference 
altitude and h c is the proper sink rate for the given glide slope. 

The digital equivalent equations are as follows: 


£ j 


Y i ~ Y ref 


0 . 

3 


k 4 " h r ) + k 5 M 9 


S = Cl + ° 2 £j “ 1 ’ 

) =e + t e h + c, e n 
£ j £ j-l 6 p j 3 p j-l 


c j 


0 . + 6 
3 ^ 
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+ 



L 



where 


c-, = e 


a£T 

t 2 


C 2 = ^ ~ C 1> • 


c 3 k 6 ^ a 3 T ' 1 ^ 


5. Flare Mode 

The flare mode controls the traditional exponential flare [25]. 
The flare has three boundary conditions, 

hp initiation altitude of feedback control law, 

h f initial sink rate, 

h^ desired vertical touchdown velocity. 

A flare law that satisfies these boundary conditions is developed 
by Neuman and Foster [3], A modified version of this flare law is 
used in our simulation. The equation that satisfies the boundary 
condition is: 

-t/a 4 

h (t) =[h f - a 4 h f ]e + a 4 h f , (26) 

o o 

and a^ is calculated as 

d 4 = ) . (27) 
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The reference sink rate is the derivative of Equation (26), 


h r = - (l/a 4 ) (h f - a 4 h f ) e 

1 o 


-t/a. 


( 28 ) 


The predictive portion of the flare law (Figure 16) has two sections, 
a step command in pitch, A0p, which causes the aircraft to begin to 
rotate, and a ramp pitch command, A0^, which begins somewhat later. 
With no other disturbance, the predictive flare commands will 
generate an approximately exponential flare. Feed-back is used to 
overcome disturbance. Equation (26) is the solution of the following 
differential equation: 


h r + a 4 (h r - h f ) = 0 , 
o 


(29) 


with the boundary conditions h r = h f at t = 0, and h r = h f at h = 0. 

0 

The feed-back version of the flare law generates a corrective signal 

when Equation (29) is not fulfilled by the actual altitude, h, and 

• • 

sink rate, h, in place of h^ and h r> The corrective signal is (see 
Figure 16): 

e fb = K f [ 1 + ^] [ h + a 4 (h - h f ) ] , (30) 


which is added to the predictive pitch command. Hence, no correction 
signal is applied when the reference path is followed. 

When the flare subroutine is entered for the first time, the sink 
rate is used to calculate decision altitudes for the predictive flare 
law commands. The altitudes at which these initial calculations are 
made are somewhat above the highest value at which the flare may be 
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started. The step command altitude is proportional to the flight-path 
angle, y ^ h/v. The ramp begins at a proportionally lower altitude. 
Since the aircraft does not begin to deviate from a straight-line 
glide path instantaneously, upon receiving the pitch step command, the 
altitude for the corrective feed-back to begin is also selected pro- 
portionally lower than the step command altitude. After these 
calculations are completed, the flare mode transfers the authority 
back to the glide-slope tracking control. 

When the step command altitude, h step , is reached, the flare 
control mode takes over completely and from the sink rate, h, calcu- 
lates and executes A6 p . Then the ramp increment, A6 R , and the feed- 
back gain constant, a^, are calculated. At this point the mode 
controller is switched to its final submode. 

In the final submode the predictive ramp pitch command is added 
to the corrective feed-back flare command. The summed signals are 
transmitted as the pitch change command, e . 

Under disturbances, the feed-back term in the flare law does not 
attempt to guide along a path fixed in space, or even hold h(t) and 
0(t) at given values 1 As long as the feed-back signal of Equation 
(30) is zero no correction is made. Disturbances, therefore, tend to 
cause translations of the touchdown point rather than large maneuvers 
to meet a given touchdown point, which would often cause hard 
landings , 
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6. Calculation of the Feed-Back Controls 

To control the flight path of an airplane automatically, it would 
be desirable to control the flight-path angle, y, directly. However, 
there is no output control variable that controls y. In linearized 

_ 4 - # 

models [1], the steady velocity, V, at which the airplane flies is 
governed by the lift coefficient, which is in turn fixed by the 

4 * 

elevator angle, implying that a constant 6^ gives a fixed V. Also, 
the flight path angle, y, at any given speed is controlled by the 
thrust in the long term, implying that the ultimate result of moving 
the throttle at fixed elevator angle is a change in y without change 
in speed. But, by physical reasoning [1], we know that initial 
response to opening the throttle is a forward acceleration, and 
initial response to elevator deflection is a rotation in pitch; hence, 
the short term and long term effects of these controls are quite 
contrary. The total picture of longitudinal control is clearly far 
from simple when we represent the aircraft motion with a non-linear 
system of equations. 

To make short and long term responses agree, the aircraft is 
stabilized in the following manner. The speed of the aircraft is kept 
nearly constant throughout the landing operation and flight path is 
controlled by means of both throttle and elevator angle deflection. 

The thrust control loop maintains constant airspeed by generating a 
thrust command signal to drive the throttle servo (Figure 17). The 
thrust command signal is derived from airspeed error, horizontal and 
vertical acceleration, and the pitch command signal. These four 
signals are processed by passing through variable gains and a 
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differencing circuit so as to generate a thrust command signal of 
correct magnitude. The elevator control loop generates an elevator 
angle command signal to drive the elevator angle servo (Figure 18) so 
as to maintain air velocity constant and control the flight path. The 
elevator angle command signal is also derived from airspeed error, 
horizontal and vertical acceleration, and the pitch conmand signal, 
and then passed through variable gains. The thrust command signal and 
elevator angle command signal are given by 

F T c = K T1 V a ' K T2 f + K T3 f + K T4 e C > 

• • 

6 E c = K D1 V a “ K D2 V " K D3 V + K D4 0 C * 

where Kj ^ and K D1 are the variable gains calculated from the system 
Equations (13), (14), and (15) under zero wind conditions, given by 


’ K n 


G 10 

G 11 

” G 12 


H 9 

K T2 
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H 10 
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H 11 
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where 


'10 


’ll 

'12 


'10 


'll 


= - (C-j + C 2 cc) cos(y' - y) - C 3 sin(Y' - y) 

= - C 4 s i n (y ' - Y) 

= Dg cos(6j + a) 

= cos(y' - y) - (C-j + a) sin(Y' - y) 

= C 4 cos(y' - y) 

= Dg s i n ( 6-p + a) 

C 5 

= C 6 

= C 7 V a cos(y' - y) + (Cg V a + Cg q + C-j g a') sin(Y' - y) 
= (Cg v a + c g q + C 10 <i‘ ) cos(Y' - Y) - Cy V a sin(Y' - Y) 
= V a +C 12 q + C 13 a' 

= GgG, - GyGg 
= G 5 G 1 - GyGg 

= G 5 Gg - GgGg 

■»f». 

• v H s 

= (G^G-j - G 4 Gg)/H 5 

* (g 2 g 7 -g 4 g 5 )/h 5 

= V H 4 

= (G 2 Gg - ^3^g)/^ 4 

= (g 2 g 6 -g 3 g 5 )/h 4 , 


with the C n values provided in the appendix. 
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CHAPTER V 


RESULTS AND DISCUSSION 

The two-dimensional equations of motion, Equations (1), (2), and 
(3), for the aircraft discussed in Chapter II were solved on a digital 
computer using a fourth-order Runge-Kutta technique. Wind models 
incorporated into the governing equations include (1) atmospheric flow 
over simulated buildings, (2) atmospheric flow in the absence of 
buildings, and (3) atmospheric flow associated with thunderstorm gust 
fronts. The influence of these wind fields on the aircraft landing 
under different conditions of terrain roughness is investigated. The 
aircraft characteristics used in the simulation are that of the DC-8 
and the DHC-6, specifications of which are given in the Appendix. The 
initialization conditions for the simulated landings of the aircraft 
with fixed controls are trimmed conditions on a -2.7 degree glide 
slope, with the descent beginning at an altitude of 91.4 meters (300 feet). 
This corresponds to a touchdown point of 1939 meters (6361 feet), 
down range from where the descent begins. Any variation in 
winds will cause the aircraft to deviate from the glide slope. The 
deviation in touchdown point for fixed control conditions is defined as 
the distance between the actual touchdown point and the intended glide 
slope touchdown point. 

Figures 19 and 20 show the descent trajectories of the aircraft 
into a wind blowing over a two-dimensional bluff-type body, similar to 
a block building, and a step, similar to a long, wide building. 
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Figure 19 Fixed control landing over a block building 



Altitude 



Figure 20 Fixed control landing over a step 






respectively. The heights of the simulated block building and the 
step building are of 20 m and 10 m, respectively. Figure 21 shows the 
trajectories of the aircraft in identical wind conditions without the 
presence of the building (i.e., the neutral atmospheric boundary 
layer), and Figure 22 shows the descent trajectories of the aircraft 
into the wind fields characteristic of a thunderstorm gust front. 

Figures 23 and 24 show the winds, both horizontal and vertical, 
that were encountered during the descent of the aircraft through the 
building-disturbed winds — a block building and a long, wide building, 
respectively. Figure 25 shows the horizontal and vertical winds that 
were encountered by the aircraft during descent through the thunder- 
storm gust front. Three flow conditions were used in each simulation 
of the wind fields for the flow over the block building, for the 
atmospheric boundary layer flow without the building present, and for 
the thunderstorm gust front. The surface roughness parameter, Z , was 
parametrically assigned the values of 0.2, 0.4 and 0.8 m, with corre- 
sponding friction velocities, u*, of 1.25 m/sec, 1.4 m/sec, and 
1.6 m/sec, respectively. These combinations of friction velocity and 
surface roughness give wind speeds of 12.3 m/sec, 11.4 m/sec, and 
10.4 m/sec, respectively, at an altitude of 10 meters. For simu- 
lation of flow over a long, wide building, surface roughness values of 
0.45 meters and 1.0 meters were used, with an assigned wind speed of 
10 m/sec at an altitude of 10 meters. This corresponds to friction 
velocities of 1.27 m/sec and 2.5 m/ser, respectively. 

Figures 26, 27, 28 and 29 show the landing trajectories of the 
aircraft with automatic controls, through the atmospheric boundary 
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Figure 22 Fixed control landing in thunderstorm gust front 
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Figure 25 Wind distribution in the thunderstorm gust front along 
the flight path 
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Figure 27 Automatic landing over block building 
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Figure 29 Automatic landing in the flow associated with thunderstorm gust front 




layer without the presence of the building, flow over a block building, 
flow over a long, wide building, and a thunderstorm gust front, 
respectively. The same wind field simulations that were used in the 
landing simulation with the fixed controls are applied here also. 
Figures 30, 31 and 32 show the controls, thrust and elevator angle, 
that were actuated by the automatic control system to track the glide 
slope through the block building, the Tong, wide building, and the 
thunderstorm gust front disturbed wind fields, respectively . Addi- 
tionally, Figure 33 shows the trajectory of an aircraft with the 
characterise cs of a DHC-6 Twin Otter with automatic controls landing 
through the atmospheric flow over a block geometry building. 

The deviations from the touchdown point in the different wind 
fields are presented in Table 1. Note that a positive deviation 
indicates a long landing and a negative deviation indicates a short 
landing. The deviation point for the automatically-controlled air- 
craft was taken as the difference between the actual touchdown point 
and the touchdown point determined by the prescribed trajectory. 

This value is computed, as illustrated in Figure 34, by adding X = 3h r , 
the specified condition for capture to begin, AX = 0.02h r required 
for the glide slope to be captured, cot y , the horizontal distance 
covered while on the glide slope, and 0.2h^ cot y^, the horizontal 
distance traversed during flare. The value of is specified as 
1 .35 degrees . 

One observes that the aircraft made short landings in almost all 
cases with fixed controls. By comparing the landing of the aircraft 
in the atmospheric boundary layer, with and without the presence of the 
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Figure 30 Controls required for landing over a block 
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Figure 31 Controls required for landing over a step 
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Figure 32 Controls required for automatic landing in thunder- 
storm gust front 
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Figure 33 DHC-6 landing through atmospheric flow over a block building 


Table 1. Deviation from touchdown point 



Wind Condition 

With Fixed 
Controls 
(meters) 

With Automatic 
Controls 
(meters) 

I. 

Flow over a building 




(1) Block geometry 




a. 

Z o =0.2m 

u*=l .25m/sec 

-198 

+35 


b. 

Z Q =0.4m 

u*=l .4 m/sec 

-309 

-30 


c . 

Z o =0.8m 

u*=l .6 m/sec 

-415 

-5 


(2) Long step geometry 




a. 

Z Q =0.45m u*=1.27m/sec 

-147 

-15 


b. 

Z Q =1 .Om 

u*=2,5 m/sec 

-185 

+11 

II. 

Flow without building present 




a. 

Z o =0.2m 

u*=l . 25m/sec 

-313 

-14 


b. 

Z Q =0.4m 

u*=l . 4 m/sec 

-328 

+7 


c. 

Z o =0.8m 

u*=l .6 m/sec 

-350 

+6 

III. 

Flow associated 

with thunderstorm gust fronts 


a. 

Z o =0.2m 

u*=l . 25m/sec 

-223 

+19 


b. 

Z Q =0.4m 

u*=l .4 m/sec 

-80 

+8 


c. 

Z Q =0.8m 

u*=l .6 m/sec 

+505 

-27 
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Figure 34 Automatic landing reference trajectory 



building, one can observe the effect of complex wind patterns caused 
by the presence of the large bluff objects. The aircraft lands 
approximately 300 meters short of the touchdown point in the atmospheric 
flow without the building, which is a direct result of the decreasing 
head wind with elevation. However, this is a predictable effect and 
all landings are within 37 meters of one another, as shown by Figure 21, 
page 55. Landing through the same atmospheric flow over buildings, 
however, causes larger variations in touchdown point between the 
different wind conditions. With a surface roughness value of 20 cm 
the aircraft lands approximately 200 meters short, whereas for a Z Q of 
80 cm the aircraft lands approximately 400 meters short. 

This variation in touchdown points can be explained by looking at 
the winds encountered by the aircraft along the flight path; see 
Figure 23, page 57, and Figure 24, page 58. The aircraft encounters 
an increased head wind (positive horizontal shear) and a downdraft 
just before the building. An increasing head wind during approach 
with fixed controls causes the aircraft to be high on the glide slope. 
This effect is absent in the undisturbed atmospheric boundary layer 
because of the continuously decreasing head wind, which causes the 
aircraft with fixed controls to always be below the glide slope. The 
aircraft does have a slight downdraft, but this is not sufficient to 
overcome the head wind influence. Just past the building, there is a 
sudden drop in horizontal wind and a sudden increase in updraft. This 
reverses the previous effect and forces the aircraft to go below the 
glide slope. Thus, competing wind effects occur relative to the 
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undisturbed wind field case, resulting in an unpredictable deviation 
from touchdown. 

One observes that the magnitude of the positive horizontal 

(tail wind) shear is largest for the case of Z q = 20 cm and smallest 

for the case of Z Q = 80 cm. Thus, landing under the conditions of 

Z = 20 cm and Z = 40 cm creates sufficiently high positive hori- 
o o 

zontal shears and updrafts to force the aircraft higher during part of 

the approach than the case with no building present, resulting in 

shorter deviation from touchdown. The case Z = 80 cm does not create 

o 

as large a positive horizontal shear and updraft and the influence of 
the decreasing head wind and downdraft in the wake of the building 
produce a shorter touchdown than occurs with the undisturbed wind 
field. 

In the case of landings through the flow over the long building, 
the aircraft lands 147 m short for the case Z Q = 45 cm and 185 meters 
short in the case Z Q = 100 cm. One should notice that in this case the 
descent starts from an altitude of about 33 meters. Therefore, these 
touchdown points cannot be compared with the touchdowns in the 
landings through the atmospheric boundary layer without the building. 

However, by looking at the winds encountered by the aircraft 
along the flight path, one can observe the similarity with the block 
building case. Just as the aircraft passes over the edge of the 
building, there is a sudden drop in horizontal wind which causes the 
aircraft to go below the glide slope, and thereby resulting in an 
unpredictable deviation from touchdown. One can observe that the dif- 
ference between the case of Z Q = 100 cm and Z Q = 45 cm is that the 
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initial head wind is higher and decreases more rapidly. Also, the 
aircraft encounters a smaller updraft than in the former case which 
tends to decrease the sink rate. Both of these wind effects result in 
a shorter touchdown, as indicated by the results in Table T, page 69. 

In the case of landings with automatic controls, the deviations 
from the touchdown points are very small for landings through the 
atmospheric flow without the building, as well as with buildings 
present. Even though the deviations from the touchdown point are not 
significant, one can see that complex wind patterns created by the 
presence of large buildings do create larger deviations in touchdown. 
Moreover, a significant factor on the safety of aircraft operations is 
the rate at which the controls must be operated to track the glide 
path. Figure 30, page 65, shows that the thrust must be reduced 
significantly as the aircraft passes over the building to compensate 
for the sudden excursions in the horizontal and vertical winds. 

Thrust was cut almost 50% in the case for Z Q = 80 cm. Figure 31, 
page 66, shows the controls that were used in the case of flight over 
the long, wide building. In this case the elevator angle changes 
quite rapidly to compensate for the effects of changing winds. 

Figure 22, page 56, and Figure 29, page 63, show the aircraft 
landing trajectories in the flow associated with a thunderstorm gust 
front with fixed controls and with automatic controls, respectively. 

In the case of fixed control landings, the deviations between touch- 
down points are large, as much as 725 meters variation from the case 

for Z = 20 cm and for Z = 80 cm. Figure 22 shows that there are 

o o J 
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very large deviations from the flight path itself and thus would 
result in unsatisfactory landings. 

Figure 25, page 59, shows the winds encountered by the aircraft 
along the trajectory. Initially, there is high vertical updraft which 
increases to approximately 13 m/sec and then suddenly drops to almost 
zero. Thus, the aircraft is initially blown above the glide slope and 
then suddenly this wind dies out, resulting in a high sink rate. In 
the case for Z q = 20 cm, the head wind remains almost constant and 
then starts decreasing very rapidly; at the same time the vertical 
updraft suddenly drops to zero, creating a strong turning moment on 
the aircraft which results in a shorter landing. For landing under 
the conditions of Z Q = 80 cm, head wind and vertical updraft remain 
constant for a longer period of time, thus keeping the aircraft above 
the glide path for a longer period of time. 

In the case of landings with automatic controls through this wind 
field, the deviations from the flight path are small, but the 
influence of these severe winds is to cause rapid changes in the 
controls. Figure 32, page 67, shows the controls that were applied to 
maintain the flight path. Here thrust starts out initially at a high 
magnitude, increases slightly and then drops very suddenly to almost 
one-third of the initial value. This is because the pitch angle is 
maintained constant in automatic landings, a vertical updraft causes a 
larger angle of attack and increased drag. Therefore, when the 
vertical updraft suddenly goes to zero, the angle of attack decreases 
and drag drops, hence thrust is reduced to prevent going high on the 
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glide slope. The elevator angle also changes very rapidly from 0 to 
-6 degrees in the process to compensate for the sudden drop in 
vertical wind and counterclockwise moment produced by the vertical 
wind, since a negative elevator setting corresponds to a counter- 
clockwise rotation with the airplane traveling in positive 
x-di recti on. 

Figure 33, page 68, shows that the aircraft with the character- 
istics of a DHC-6 on a 6-degree glide path behaves basically the same 
as a DC-8. For the landing of this aircraft, few changes were 
required in the automatic control system. In the glide-slope capture 
mode the reference flight path angle was changed from 2.7 to 6 degrees 
and the gain constant, k^, from 0.05 to 0.1. In the glide-slope 
tracking mode the constant, a 3 , is changed from 0.01 to 0.02, and in 
the flare mode the constant, a^, is changed from 4.0 to 5.0. These 
changes in different control modes increase the response of pitch 
angle command variable, 0 , 
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CHAPTER VI 


CONCLUSIONS 

The two-dimensional aircraft landing simulation study has pro- 
vided some basic results concerning the problems of wind shear due to 
the presence of large buildings or other bluff geometries. The air- 
craft encountering a strong wind shear caused by the edge of a 
building is drawn towards the building. With fixed controls, 
deviation in touchdown point in excess of 400 m resulting from 
variation of the horizontal wind during the final 100 m of descent has 
been computed under wind shear conditions that may realistically be 
encountered around buildings or bluff terrain features. 

Wind shear due to thunderstorm gust fronts can cause very severe 
departures from the glide slope during landing. Although thunderstorm 
gust fronts are not encountered too frequently, landing through such a 
gust front can be very hazardous and requires rapid changes in the 
controls required to maintain the flight path. Based on the model of 
thunderstorm field investigated in this report, changes in the controls 
at the rate of 7 degrees of elevator angle and 9072 kg of thrust per 
one-half second for the DC-8 were required. 

Although the surface roughness parameter, Z , shows little 
influence on touchdown points during the landings through the atmos- 
pheric flow over level terrains, it does have considerable effect on 
aircraft landing when large buildings or bluff objects are present 
near the airports. 
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APPENDIX 


The aircraft characteristics used in this simulation study are 
for that of a DC-8 and a DHC-6. The initial conditions and aircraft 
physical data are given as follows [7]: 


H(m) 

V a ( m/sec) 

y(deg) 

W(kg) 

I yy (kg -"’ 2) 

L-j-(m) 

S T (deg) 

c(m) 

s(m 2 ) 


DESCRIPTION DC-8 

Reference altitude 91 

Initial velocity 70 

Initial flight path angle -2.7 

Aircraft weight 90700 

c 

Moment of inertia 5.3 x 10 

Moment arm of thrust vector 1.2 

Angle between F T and FRL 3.15 

Chord length 7 

Wing area 256 


DHC-6 

91 

46 

- 6.0 

4985 

3.2 x 10 4 
-0.91 
0.0 
2 
39 


The expressions for the aerodynamic coefficients of the DC-8 
aircraft are: 


C, = 


Cn = 


C. + C, a' + C. <5 r + -59- C. + — C. , 
L- L - L 5 E 2 V* L„ 2V. L- ’ 


o a 


^D + ^D a ' + 2 a ' * 

o a a 


E a q L, a a 

2 


C = 
m 


C + C a' + C c + c 

m. m fi E 2V a m n 2V a m-, ’ 

°E 


o a 


a q a' 
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where the usual notation is used for the various stability derivatives. 
The values of the stability derivatives are given as follows: 


DC-8 DHC-6 



0.90 

0.86 

C L 

a 

5. 30/ rad 

6.109/rad 

C L 

6 E 

0.0053/deg 

0.5236/deg 


7.68/rad 

2 . 1 52/rad 

C L. 

a 

0.0 

0.0 

C °o 

0.140 

0.32 

C D 

a 

0.501/rad 

0.9832/rad 

C D 2 

a 

1 .818/rad 2 

0.0 

o 

E 

CJ 

-1.01 

0.0/rad 2 

C m 

a 

-1.062/rad 

-2.026/rad 

n 

L/ 

rn ~ 
S E 

-0.0161/deg 

-2.068/deg 

m q 

-12.30/rad 

-28.76/rad 

m* 

-4.01/rad 

-8.663/rad 


a 


The various dimensionless groups 
study are provided below: 

D, = 6 - — 

1 2 w 


D = 

Up o 

V a 


(normalizing factors) used in this 
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D = J£ 

3 v 2 

a 


D = — 
u 4 2H 

n - ° SCH 2 
°5 ' 2 g I 


D 


= JitL 


6 ,,2 
w V a 


L x Do 

n = 

U 7 I 

yy 


The following are the "C" coefficents used in variable gain computations of 


F t and 6 r : 

l r t-p 


c i D 1 C D 


C 2 D 1 C D 2 

a 


C 3 D 1 C L 


C 4 = D 1 C L 


Cr = Dr C 

5 5 m 


Cr = Dr C 

6 5 m. 


C 7 °1 C D, 


C 8 D 1 C L, 


C 9 " D 1 D 4 C, 
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C 10 = D 1 C L* 
C n = D 5 C m 0 
C 12 = °5 C|n q 
C 13 = °5 °4 C m-. 
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